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Abstract 

A robust implementation of a Dupire type local volatility model is an 
important issue for every option trading floor. Typically, this (inverse) 
problem is solved in a two step procedure : (i) a smooth parametrization 
of the implied volatility surface; (ii) computation of the local volatility 
based on the resulting call price surface. Point (i), and in particular how 
to extrapolate the implied volatility in extreme strike regimes not seen 
in the market, has been the subject of numerous articles, starting with 
Lee (Math. Finance, 2004). In the present paper we give direct analytic 
insights into the asymptotic behavior of local volatility at extreme strikes. 



1 A new formula for local volatility extrapola- 
tion 

Volatility remains a key concept in modern quantitative finance. In particu- 
lar, the Black-Scholes implied volatility surface ctbs = '^bsI^^tT) is the central 
object of any option trading desk, see e.g. [TS]. On the quantitative and com- 
putational side, a smooth and arbitrage free parametrization of the implied 
volatility surface is a crucial step towards a robust implementation of a Dupire 
type local volatility model [HI [12]. Indeed, Dupire's formula 



implies that any arbitrage free call price surface 

C = C{K,T) = CBs{K,T;aBs{K,T)) 

which arises from a (not necessarily Markovian) Ito diffusion is obtained from 
the one-factor ("Dupire's local vol") model 

dSt/St = <Jioc{St,t)dWt. 



partially supported by MATHEON 



(Note that spot remains fixed in the present discussion, and that we work under 
the appropriate forward measure to avoid drift terms.) It is helpful to think 
of local volatility as a Markovian projection, term coined in [22], of a higher 
dimensional model (e.g. Heston); the first component then forms an Ito diffusion 
of the form 

dSt/St = astoch{t,i^)dWt. 

Indeed, it is known (see e.g. [H] and the references therein) that a^^^{K,T) = 
E[(7gjQ^jj \St = K\\ in practice, this means that even for stochastic volatility 
models with fully explicit Markovian specification, sampling from the corre- 
sponding local volatility models requires substantial computational effort. (In 
particular, the singular conditioning requires Malliavin calculus techniques, as 
was pointed out e.g. in [T7].) 

The analysis of implied, local, and stochastic volatility and their interplay 
has been subject of countless works; a very small selection relevant to the present 



discussion is [TJ [H 31 HH [121 [20] ■ Our contribution here is a formula (( 1.6 ) below) 
that allows for approximation of a'^^^{K^T) when K is large (and similarly, K 
is small). The main ingredient to this formula is a known moment generating 
function (mgf) of the log-price {Xt) (under the pricing measure), 

M{s,T) exp(m(s,r)) :=Eexp(sXT), 

assumed to be finite in some (maximal) interval (s_(T), s_|_(T)) with critical 
exponents s_ and s+ defined as 

s_ (T) := inf {s : M(s, T) < oo} , s+ (T) sup {s : M (s, T) < oo} . 

We also assume that call prices have sufficient regularity to make Dupire's for- 



mula ( 1.1 1 well-defined, and that the mgf blows up at the upper critical moment: 



lim M(s,T)=cx). (1.2) 

sts+(T) 

This holds, e.g., in the Heston model p^, with log-price Xt — \og{St/ So), where 

dSt^St^tdWt, S'o = so>0, 

dVt ^{a + bVt)dt -f c^JVtdZt, Vb wo > 0, 

with a > 0, & < 0, c > 0, and d {W, Z)^ = pdi with p e (-1, 1). We wiU prove 
the following theorem, which is reminiscent of Lee's formula |20j for implied 
volatility. 

Theorem 1. In the Heston model with p < (the relevant regime in practice), 
the following local volatility approximation holds^ 

lim '-^^ - , \ , ^ (1-3) 
fe^oo k T s+{s+-l)Ri/R2 ^ ' 



^ By a common abuse of notation, we write a^^^{k, T) instead of o'j^^^(e'"", T) when we wish 
to express the local vol as a function of log-strike k. 
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(1.4) 



where k — log (K/Sq) , s+ — s^(T) and 

= c^s+{s+ ~ 1) [c2(2s+ - 1) - 2pc{s+pc + b)] 
- 2{s+pc + b) [c^{2s+ - 1) - 2pc{s+pc + b)] 
+ Ape [c^s+{s+ - 1) - {s+pc + bf] , 
i?2 - 2c\s+{s+ - 1) [c2s+(.s+ - 1) - [s+pc + . (1.5) 

The origin of this resuh hes in the saddle point based approximation formula 



2^m{s,T) 



s(s-l) 



(1.6) 



s{k,T) 



where k denotes log-strike, and s — s{k,T) is determined as solution of the 
equation 

^m[s,T)=k. (1.7) 



As such, our formula ( 1.6 ) is not restricted to the Heston model. As a trivial ex- 



ample, let us consider the generalized Black-Scholes model with time-dependent 
volatility, 

dSt ^ St^/W)dWt. 
We find m{s,T) = \s{s — 1) v{t)dt, and then, correctly, 

aL(fc,T)=«(T). 

(In this example the evaluation of s{k,T) plays no role, since the fraction on 
the right hand side of (1.6) does not depend on it.) 



In fact, we expect our approximation formula (1.6) to work whenever the 



saddle point method is applicable (also assuming that call prices are smooth 



enough to make (1.1 1 a well-defined quantity); the essence of the argument is 
given m Section ^~0i course, the ultimate justification of a saddle point ap- 
proximation involves tail estimates which may present mathematical challenges 
(while easy to observe numerically); in the Heston case, we achieve this by a 
subtle application of ODE comparison results, applied to the underlying Riccati 
equations (cf . Appendix |A]) , thus completing our proof of the above theorem. 
The asymptotic equivalence of (1.3) and (1.6) is discussed in Section [4] 

Interestingly, even when the blow-up of the mgf is too slow to apply the 



saddle point method, the approximation formula (1.6) can give surprisingly ac 



curate results. Our attempt to understand this phenomenon, besides numerical 
evidence in the variance gamma model for suitable parameters, discussed in 
Section (l5), passes through Karamata's Tauberian theorem and is the content 
of Section 3j We have not pushed our investigations too far, however, since the 
meaning of Dupire's local volatility in the presence of jumps may be questioned 
(cf. our comment below on extension of Dupire's formula to jump settings.) 
Various additional comments are in order. 
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1. Equation (1.7) is solvable for large k, since (1.2) implies 



lim —m(s,T) — oo. 

sts+ OS 



2. If there is no blow-up, i.e. (1.2) does not hold, then (1.6) is typically 
incorrect. See Example |4] in Section [5] for some hints on how to handle 
such cases. 

3. We have s{k,T) f s+{T) as fc — > oo; hence, in models with moment 
explosion [Tl I20j. where s+(T) < oo, the denominator in (1.6) may be re- 
placed by s+{T){s+{T) — 1). While this is correct to first order, it is often 



6. 



preferable to use (1.6) as it is, and to calculate s(fc,T) by (numerically) 



solving (1.7). This tradeoff between simple formulas and numerical preci- 



sion is illustrated in several examples in Sections [4] and [Sj The comment 
applies in particular to the Heston model. 



4. There is a version of our approximation formula (1.6) for small values of 
K (i.e. K I 0, 01 k I — oo), which requires that the mgf blows up at the 



lower critical moment s_(T). If fc < and |fc| is large, equation (1.7) has 
a unique solution s_(fc,r) < 0. Then the approximation (1.6) holds, if s 
is replaced by s_. 

5. There are extensions of Dupire's work to jump diffusions and also pure 
jump models; the resulting "local" version of these models is studied in [3]. 
Local Levy models were introduced earlier in [5]. In particular, Dupire's 
formula (which may be written as a PDE) becomes a PIDE which features 
an integral term involving the second derivative of C w.r.t. strike, times 
a kernel depending on K, integrated against all strikes in (0, oo). (The 
formula, which we need not reproduce here in full technical detail, appears 
in Theorem 1 of [3].) 

Another difficulty in the jump setting is the potential lack of immediate 
smoothing. For instance, the variance gamma model satisfies the above 
PIDE only in viscosity sense; in fact, call prices in the variance gamma 
model may not be twice differentiable in K for small times, as was noted 
in But for sufficiently large times our formula (1.6) works, see Exam- 
ple [3] in Section [5j 

We conclude that, in a general jump setting, Dupire's formula, as stated 



in ( 1.1 ), may be ill-defined; moreover, even if call prices are smooth enough 



to make the formula well-defined, it fails to recreate the correct marginals 
of the original price process. 

Even so, given the industry practice of applying Dupire's formula to any 
given call price surface, we discuss in Section [5] what happens when ap- 
plying (1.6) to jump models, if possible. Formula (|1.6|) simplifies in expo- 



nential Levy models, which have the property that m(s,T) is linear in T; 
thus, the numerator in ( 1.6 ) may be replaced by 2m(s, 1). In jump models. 
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we also expect crf^ci^,T) to explode as T | 0, and we shall observe and 
quantify this blow-up in some examples below. There is potential practi- 
cal value in that a Dupire local volatility surface, fitted to market data, 
may so be inspected for evidence of jump behavior (thereby questioning 
the use of Dupire's formula in the first place). 



2 Saddle point asymptotics 

As is well known [7], we can recover the call price C and the probability den- 
sity D{-,T) of St by Laplace-Fourier inversion from the mgf: 



e 



k 



-ksM{s,T) 



C{K,T) = — e-'^-^ds, (2.1) 
^^^^ J-ioo s{s-l) 



too 



D{x,T) = —l e^(''+i)'°s^M(s,T)ds. (2.2) 



— 'too 



2iTT 

Now differentiate the call price under the integral sign: 

drCiK, T)^^ r ^^T^e-'^Uf (., T)ds. (2.3) 
2«7r J_,^ s(s - 1) 

By Dupire's formula, we have 

' ^ " K^D{K, T) - p_^^ ,-ksMis, T)ds ' ^ ^ ^ 



Both integrands in (2.4 1 have a singularity at s = s+, since M{s,T) gets infinite 
there. The singular behavior of AI{s,T) dominates the asymptotics of both 
integrals. The resulting asymptotic factor cancels, and only the contribution of 



2 3rrn,(a,T) pgjj^g^j-^g Tj^jg jg t]^g jjea behind (|1.6|. 
s(s— 1) 1 y 



s(s- 

To implement it, we analyze both integrals in ( 2.4 ) by a saddle point approxi- 
mation [TU] . If M features an exponential blow-up at the critical moment s+ , its 
validity can be justified rather universally. Examples include the Heston model, 
double exponential Levy, and Black-Scholes. (Note that the critical moment is 
s+ = oo for Black-Scholes.) If the saddle point method is not applicable (be- 
cause of insufficient blow-up), different arguments are required; see the following 
section. 



So let us proceed with the saddle point analysis of (2.4). For both integrals, 
we only use the factor e~'^*M(s) to find the location of the (approximate) saddle 
point s = s(fc,T). The saddle point equation is ( |1.7[ ), obtained by equating the 
derivative of e^^''M{s) to zero. We move the integration contour through the 
saddle point. Then, for large fc, only a small part |9(s)| < h{k) of the contour, 
around the saddle point, matters asymptotically. (The choice of the function h 
depends on the singular expansion of Af; see Section |4] for an example.) The 
integral can be approximated via a local expansion of the integrand. Let us carry 
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this out for the denominator of ( |2.4[ ). (In the following formulas we write m" 
for d'^m/ds^.) 



s+ioo 



e-^'M{s, T)ds 



s — ih{k) 



^-ki 



M{s,T)ds 



S'\-ih{k) 



ih{k) 



exp {^ks + m(s, T) + k{s — s) -\- ^ra" {s^ ^)(^ ~ ^)'^) ds 



^m{s,T) — ks 



s-\-ih{k) 



exp{^m"{s,T){s- s)^)ds. 



(2.5) 



' s — ih{k) 

In the Taylor expansion of the exponent we have used the equation m'(s, T) — k. 
Now the crucial observation is that the numerator of (2.4) admits a similar 

?Tm(s,T 



approximation, where the only new ingredient is the factor 2 ^T'm.{s,T) _ 



<s ~ 1) 



e'''"'M{s,T)ds 



s~ih(k) s{s — 1) 

~ 2e'"(^^^)-'=^ ^^^^ (1 + 0(1)) exp (im"(s, r)(s - s)^) ds 

J~s-ih{k) - 1) 

1 /' ^ T^A ps-{-ih(k) 

^ 2^^P^e™(^'^)-'^- / exp(im"(s,T)(s-s)2)ds. (2.6) 

Dividing (2.6) by (2.5) concludes the derivation. Summarizing, we note that 
the asymptotics of o^(fc) are governed by the local expansions at s = s of the 
integrands in (2.4). The respective first terms of both expansions agree, and 
thus cancel, except for the factor (1.6). 



3 Algebraic singularities and Karamata's theo- 
rem 

The saddle point method is well suited to treat mgfs of exponential growth, 
such as M{s, T) k, exp(l/(s+ — s)), but fails in cases of slower blow-up. To see 
how to analyze these, let us assume that the mgf M grows like a power at the 
(finite) critical moment: 

(s+ - sY^ 

(The variance gamma model is a typical instance.) The quantities ci = ci{T) > 
and C2 = C2{T) > are independent of s, but may be functions of matu- 
rity T. (In particular, we assume that C2 does depend on T, which holds in 
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Levy models.) Since 



the saddle point s satisfies 



ds 



m{s,T) 



C2 



k ■ 



Inserting this into the time derivative 

d 



TO(s,T)^C2(r)log — 



of m yields 



m (s, T) 



s{s-l) 



s=s{k,T) 



2c2(r) logfc 

s+{s+ - 1) ' 



(3.1) 



To justify ( 1.6 ), we now have to argue that (jf^^{k, T) has the same asymptotics. 
Again, we use the representation (2.4). To put it briefly, the reason why the 
approach from the preceding section fails is that one cannot find a suitable h{k). 
(Either the tails 19(5)1 > h{k) of the integrals are not negligible, or the local 
expansion is not uniformly valid.) But it is still true that the local behavior of 
the integrands near s+ fully determines the asymptotics of the integrals in (2.4). 

We write f{-,T) for the probability density of the log- price Xt- Note that 
the denominator in (2.4) equals 2i7r/(fc,T). The (one-sided!) Laplace transform 
of fc f{k,T) satisfies 



e-'''e'+''f{k,T)dk^cis''^, 



This follows from 



e'^f{k,T)dk+ / e'''f{k,T)dk. 
-oo Ja 



siO. 



Cl 

{s+ - 



(3.2) 



(3.3) 

since the first integral in (3.3) is 0(1). By Karamata's Tauberian theorem (Sj 
Theorem 1.7.1], we obtain from (3.2) that 



'f{x,T)dx 



cik" 



r(c2 + i)^ 



k^ 



hence, by differentiating. 



fcCl 



r(c2) 



Similarly, the asymptotics 
dTm{s, T) 



s{s-l) 



-M(s,T) 



C2 



log 



OO. 



(3.4) 



sts+, 
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imply that the numerator m (2.4) approximately equals 



e-+'=fc--i log k. (3.5) 

s+(s+ - i)r(c2) 



Now divide (3.5) by ( |3.4[ ) to see that o'j^Qj.(fc,T) approximately equals (3.1). 
Note that we did not talk about Tauberian conditions, which are necessary to 
make this derivation rigorous, such as monotonicity of the density. In concrete 
cases, where an analytic continuation of the mgf is available, a Hankel contour 
approach [T3] might be preferable to Karamata's theorem. 



Local vol at extreme strikes in the Heston 
model 



In this section we give a numerical example and explain how (1.3) is obtained 



by specializing (1.6). (But recall that (1.6) is so far just a recipe and not a 
theorem; a rigorous proof of (1.3) is given in Appendix [A|) 

Figure[l]compares the approximations ( 1.3 ) and ( 1.6 ) for the local vol. While 
asymptotically equivalent, the plot suggests that (1.3) has an 0(1) error term, 
whereas the error of (1.6) seems to be only o(l). Note that the right hand side 



of ( 1.6 ) can be easily evaluated numerically, by using the explicit expression 



of the Heston mgf in ( 1.7 ). We will now show that the right hand side of ( 1.6 ) is 




Figure 1: Local volatility squared o'^Q^(fc,r) (solid curve) in the Heston model, 
with parameters T = 1, a = 0.0428937, b = -0.6067, c = 0.2928, Sq = 1, 



vo = 0.0654, p = —0.7571. The approximation (1.6) is dashed, and (1.3) is 
dotted. 



indeed asymptotically equivalent to the right hand side of (1.3). This requires 
us to know that 

^ ' - ^2k/a, (4.1) 

s=s(k,T) 



2 QfHs,T) 
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where a = (t(T) is the so-called critical slope, defined as 

dT* 



a{T) 



ds 



(4.2) 



r*(s) sup{t > : £[6"-^'] < oo} 



In fact, while the computation of the critical exponent s_|_ in the Heston model 
requires simple numerics, the critical slope can be computed in closed form |14] : 
we have a{T) = TR1/R2, where Ri — Ri{b,c, p, s+(T)), i — 1,2, are defined 
in ([L4l)-([L5l). 

Since s{k, T) — > s-|_(T) as fc 00, the right hand side of ( 1.6 ) then satisfied 



2^m(s,T) 



s(s-l) 



X fc, 



00, 



^iT)s,iT)is,iT)-l) 

which is the formula from Theorem [T] Let us now discuss validity of (4.1 1. The 
argument which follows nicely illustrates how formula (1.6) is used in stochastic 
volatility models of affine type. First, m{s,t) « VQip{s,t) for a function -0 for 
which we knovijf] 



and also 



— K 1 



ttT*{s). 



If we write s+ = s+(r) when T is fixed, this translates to 

""0 . 



d_ 

d_ 
df 



m(s, t) 
m{s, t) 
m{s, T) 



'+1 



Equation (1.7) leads to s 



\(j{s+ - s) 

fCT(s+-s)2' 

. , ^t^- 

fK.+ -,s))2' 

pk-^/"^ + o{k-^/'^), since 



(4.3) 
(4.4) 



d_ 
d~s 



m{s, t) 



with /3 — Substitution then yields 

d 



/3k- 



^1/2 



dT 



m{s,T)U 



-a^f/k 



k/a, 



which concludes our derivation of (4.1). 



■^It is worth noting that cr(T) ~ const X T as T — > 00. This suggests that a'^^^{kT,T) 
admits a non-degenerate hmit as T — > 00; Gatheral's SVI limit of Heston implied volatility 
was obtained in a similar regime. 

•^This follows from a straightforward analysis of the Riccati equations |14j . 
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5 Some remarks on Dupire's formula for jump 
models 

As discussed in the introduction, a direct application of Dupire's formula is not 
easy to justify in the presence of jumps. Even so, given the industry practice of 
applying Dupire's formula to any given call price surface, we now discuss what 



happens when applying formula (1.61 to some examples of jump models 



Example 2 (Double exponential Levy). For zero drift, the mgj is given by JEI 

x+p , x-ii -pY 



M(s,r) =exp T 



+ A 



A. 



The critical moment is s+ = A.)-, and the saddle point is located at 



XX+pT 



Formula (1.6 1 thus yields 



2y/Xp 



v/A^(A+ - 1) 



(5.1) 



(5.2) 



In Figure^ the fit of (5.2| is not satisfactory (the dotted curve) . Similarly to the 
Heston model, the approximation (5.2 1 has on 0(1) error term, whereas (1.61 
seems to have o{l), and gives a very good estimate. 




0.06 



0.04 



0.02 



Figure 2: Local volatility squared o'foc(^'^) (solid curve) in the double expo- 
nential Levy model, with parameters T = 1, a = 0.2, X = 10, p = 0.3, A_ = 25, 
A+ = 50. The approximation (|1.6|) is dashed, whereas (5.2 1 is dotted. 



As mentioned in Section^ one expects local volatility to explode for T \. 0; 
We have crf^^{k,T) « const x T^^/"^ in the double exponential Levy model. 
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To see this, note that in Levy models the saddle point s{k,T) is a function of 
k/T, and that the saddle point method works for T I as well as for fc — >■ oo. 
Therefore, (5.2) is true for fixed k andT^O, too. 



Example 3 (Variance gamma). The mgf is given by I21j 



M{s,T) 



1 



1 — 9vs — ^a'^vs'^ 



We assume that T > v /2, which guarantees that the log-price has a density, and 
hence that call prices are (see Example 1 in J3l)- The critical moment is 



Since 



T 1 

m(s, T) = log M (s, T) log , 

V s+ ~ s 



the saddle point satisfies 



By (1.6 1, we thus have 



T 
uk 



21og(fc/T) 
vs+{s+ - 1)' 



(5.3) 



According to Figure this approximation kicks in only for fairly large values 
of k. As in Section^ and Example^ an improved estimate is obtained by 
using (1.6 I directly. 



0.08 - 




Figure 3: Local volatility squared crfocC''' ^) (solid curve) in the variance gamma 
model, with parameters T = 1, ct = 0.261652, 6 = -0.218033, y = 0.0552584. 



The approximation (1.61 is dashed, and (5.3 1 is dotted 
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Example 4 (Normal inverse Gaussian). This is an example where condition ( 1.2 ) 
is violated, and our formula (1.6) does not hold. The mgf 



M{s, T) = exp (^T (^a^ - (3^ - ^ o? - (/3 + s)- 



has no blow-up at the critical moment 

s+ = a- f3, 

but a square-root type singularity, with local expansion 



M{s,T) « e*^^^^ (l _ JT^y^T^ 



(5.4) 



It is still true that o'i^^{k,T) asymptotically depends, via (2.4|, on the local 
behavior of M{s,T) near s+. However, the approximation (1.6) hinges on the 
first term of the local expansion of M{s,T). It therefore fails to capture the 
asymptotics of a^^^{k,T), which depend on the first singular term (the term 



Y^s+ — s in (5.4) ). An analysis can be done by a Hankel contour approach \13^ . 
and yields 

ST^/c 



2 1 



.{k,T) 



The numerical fit is not very good, though, and further terms should be computed 
for improved accuracy. We propose to return to this model and the more general 
GH (generalized hyperbolic) model in a future study. 

The fact that crf^^[k, T) converges to a constant might be understood by com- 
paring the NIG marginals with those of Heston's in the time T — >■ oo regime ( this 
link is made precise in fl9f ). In particular, the result is then consistent with the 
Heston asymptotics (1.3) of local vol; note that the right hand side of (1.3) is 
0(l/r) forT^oo. 



6 Conclusions 

We propose a new formula that expresses local volatility for extreme strikes as 
a computable function of commonly available model information. In the Heston 
model this leads to a proof that local volatility (squared) behaves asymptotically 
linear in log-strike (which is qualitatively similar to Lee's result ^20j for implied 
volatility). 

Although we suspect that this Lee-type behavior remains true for models 
similar enough to Heston (e.g. local stochastic volatility models with a Heston 
backbone |17j). qualitatively different behavior is seen in models with jumps. 



We derived this by applying our generic approximation formula ( 1.6 ), supported 
by numerical examples. 

While this enhances our knowledge of local volatility in a variety of models, it 
also has a clear impact on calibration of local volatility to market data. Indeed, 
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liquid option data is typically available only in a restricted range of strikes and 
maturities; our results can then be used to extrapolate local volatility in a way 
that is consistent with Heston stochastic volatility or other chosen models. In 
particular, this approach avoids any arbitrage possibilities introduced by ad- 
hoc specifications of the implied volatility surface. We also believe the present 
methodology will turn out useful in the calibration of local stochastic volatility 
models to market smiles. 

Acknowledgment. We thank Rama Cont and Jim Gathcral for helpful 
discussions. 



Proof of Theorem [l] (local vol approximation 
for the Heston model) 



By the exponential decay of the Heston mgf towards ±ioo, the formulas (2.1 ) 



(2.4) are correct for the Heston model. For the saddle point analysis of (2.4) 



we employ the approximate saddle point 

^approx(^) • — ^ 



obtained by using (4.3) in (1.7) 



(Recall that j3 = 



\/2wo/ 



cy^, and that a 

denotes the critical slope defined in (4.2).) This approximate saddle may be 
used for both integrals in (2.4). As for the denominator, this was carried out in 



detail in |14| . where an expansion of the Heston density was determined. The 
analysis of the numerator in (2.4) is similar, except that a new tail estimate is 



required. But first we discuss the local expansion around the saddle point. Let 
us fix a number a G (|, |) and define h{k) — fc^". Then, in the central range 
|s - Sapprox(^)| < K^) , we havc 



1 



1 



s(s-l) 



s+(,s+ - 1) 

1 



1) 



and (cf. (4.4)) 



2/3^ 



cr(s+ - s)2 



O 



+ 0{s+ - s) 

i + o(fc-i/2; 



1 



2/3' 

a 
2k 



(l + 0(fci/2-«)). 



Therefore, the local expansions of the two integrands in (2.4) agree, up to a 
factor that is given by 



2k 



2dTm{s,T) _ 

S(S-I) CTS+(S+-1) 



(l + 0(fci/'-")) 



(A.l) 
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where the error term holds uniformly w.r.t. the integration variable s. According 
to Theorem 1.2 of [H], we have 



1 

2iTT 



'■'M(s, T)ds ~ ^ie(i^'^^)'=+'^^^fc-3/''+'^/'=' (A.2) 



for certain constants Ai, A2 = 2/3, and = + 1. Analogously, we derive 
from (lA.ip that 



1 
2i7r 



ih{k) 



sis-1) 



e-^'M{s,T)ds 



2k 



as+{s+ - 1) 



X Aie(i-^^)'=+^^^A:-3/4+«/-'. (A.3) 



Dividing ( A.3[ ) by (A.2) shows our claim (1.3 1, provided that the tails |s — 
Sapprox(fc)| > h{k) of the integrals can be discarded. For the denominator 
of (2.4), this was shown in Lemma A.3 of Tf . So we proceed with the numera- 
tor. We consider only the upper tail, as the lower one is handled by symmetry. 
By Lemma A.3 of p3|, there is a constant B > such that 



e-'"'M{s,T)ds 



< e(i-^^)'= exp(v42%/fc- ^p-^k^/^-^^ + Oilogk)). 



(A.4) 



From (4.3) we obtain 



dTm{s, T) 



sis - 1) 



< const X k 



for all s on the contour in ( A.4[ ). This estimate can be absorbed into the factor 
exp(0(logfc)) in (A.4), so that we conclude 



< e'-^-^'^'' exp(A2 Vfc - i/3-^/c3/2-2a ^ q^^^^ ^^ 5) 

This grows slower than (A.3) (compare the relevant factors fc-3/4+a/c ]-ggp_ 
exp(-i/3"^fc^/2~2")). As for 5(s) > B, it was shown in [14 (Lemma A.2) that 



e-'"'M{s,T)ds 



0(exp((l- A3)A: + /3^)). 



This was deduced from the exponential decay of M(s, T) for large 3(s) (Lemma 
A.l in [HI). The following lemma implies that the new factor dTm{s, T)/{s{s — 
1)) grows only polynomially, so that the exponential decay of the integrand 
persists for the numerator of (2.4). This finishes the proof of Theorem [T] 
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To state the lemma, recall that m(s, t) — (f>{s, t) + vo'>jj{s, t), where (j) and ^ 
satisfy the Riccati equations 

= ai;, 0(0) = 0, 

^(s^ - s) + \(?^'^ + 6Vj + spcV', i/'(0) = 0. 

We have to show that to grows only polynomially as 9(s) — oo. Because of the 
Riccati equations, it suffices to show this for i]). Let us write ip = f + ig and 
s = £,+iy. 

Lemma 5. Let T > 0, and assume that the real part of s stays bounded in some 
interval 1 < £, < Cmax- Then, there are positive constants Ci^T fi = 1,2,3,4) 
such that for y > yo, where y^ depends only on ^nj^x '^'^'^ the other (fixed) model 
parameters of the Heston model, 

-C^.Ty" < fit) < -Ci.Ty. 
-Ca^tV^ < git) < C2,TV- 



In fact, we can take 



Ci^T = 1/ (3c) , 

C2,T 2 ^ 1) T, 

C3.T = T 1 + 



3,T — [ i + "2" 2,T 
C4,T = '2,C^tT(?C2.T- 

Proof. It follows from the proof of Lemma A.l in [Uj that (e.g. with Ci^t '■= 

fit) < -T9y = --^y <-^y=: -Ci^tV- 
c Ac 

We next provide a similar upper estimate for g. To this end we first show that 
g = git) remains > for all times < > 0. The differential equation for g, 

9=\i'2iy-y) + c^fg-ig, 5(0) = o, 

implies the first order Euler estimate 

5 (0 = 3(0) + I ^ i^iy -y) + c'f iO)giO) - 75(0) 1 1 + o{t) 

= li^^y-y)t + oit), 

" V ' 

>o 

and hence g is positive (even strictly so) on some interval (0,£i). Assume this 
interval is maximal in the sense that g(£i) = and g is (strictly) negative on 
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some further interval {61,62)- Clearly then g{6i) < 0, which contradicts the 
information from the differential equation: indeed, using (?(ei) = 0, we obtain 
the contradiction 

>o 

The observation that g > is useful to us, since it leads, together with / < 
— Ci,t2/ and 7 > 0, to the differential inequality 

9 = ^ i^^y -y) + (?f9 - 19 

<\my-y)-{c^c^,T + i)9 
<\{'^iy-y), 

and hence to the upper estimate 

VO < t < T : g{t) < \ {2^^^ -l)Txy=: C^^ry. 

We can feed this upper estimate on 9 back in the differential equation for / to 
obtain a lower estimate 

f = \[e-y''-i)^i^[f-9^)-if 

>\{e-y'-i) + "^f-"^Clry'-lf 
= -\{l + c'Clr)y' + \{e-^)-lf+'^f 

>-\{i + c'cl^)y' + \{.e-^)-if 
> - (1 + "^cl^) y' - If, 

where in the last step we assume that y is large enough so that the extra amount 
subtracted (at least: \y'^) is larger than ^ (^^ — which remains bounded. We 
also know that f{t) < -Ci,t2/ < for all < i < T. It follows that -7/ > 
and omission leads to our final lower bound on /, namely 

f>-{i + jClT) y'- 

This entails immediately 

fit) >-t(i + ^Cl^ ) y' =: -Cs,Ty'. 
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At last, wc need a lower bound on g. Again, we look for a suitable differential 
inequality. Since 5 > 0, 

9=\ (2^y -y) + c^fg - 79 

>li2^y-y)~c^\f\9-ig 

for y large enough such that the additional subtraction of Cs^tU^c^ takes care 
of 7. Using the upper estimate on g (linear in y), and the fact that (2^ — 1) > 0, 
we conclude 

g > -2C3,t2/^c^C'2,t2/- 
It immediately follows (since ^(0) = 0) that 

yO<t<T: g{t) > -2C3,TTc^C2,Ty^ =: -Ci^rV^. 
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